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ABSTRACT 



Studies of molecular outflows in high-mass young stellar objects reveal important information about the formation process of massive 
stars. We therefore selected the close-by IRAS 17233-3606 massive star-forming region to perform SiO observations with the SMA 
interferometer in the (5-4) line and with the APEX single-dish telescope in the (5^) and (8-7) transitions. In this paper, we present 
a study of one of the outflows in the region, OFl, which shows several properties similar to jets driven by low-mass protostars, such 
as HH21I and HH2I2. It is compact and collimated, and associated with extremely high velocity CO emission, and SiO emission 
at high velocities. We used a state-of-the-art shock model to constrain the pre-shock density and shock velocity of OFI. The model 
also allowed us to self-consistently estimate the mass of the OFl outflow. The shock parameters inferred by the SiO modelling 
are comparable with those found for low-mass protostars, only with higher pre-shock density values, yielding an outflow mass in 
agreement with those obtained for molecular outflows driven by early B-type young stellar objects. Our study shows that it is possible 
to model the SiO emission in high-mass star-forming regions in the same way as for shocks from low-mass young stellar objects. 

Key words, stars: formation, stars: protostars, ISM: jets and outflows, ISM: individual objects: IRAS 17233-3606 
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1. Introduction 

■ Two main theoretical scenarios, based on accretion, are usu- 
' ally proposed to explai n the formation of O - B typ e stars: 1) the 
. core accretion model (McK ee & TanL I2002L 120031) . where mas- 
sive stars form from massive cores, and 2) the competitive ac- 
cretion model (iBonnell et aU l2007h . where a molecular cloud 
fragments into low-mass cores, which form stars that compete to 
accrete mass from a common gas reservoir Both models predict 
the existence of accretion disks around massive young stellar 
objects (YSOs), and jets driving molecular outflows. The core 
' accretion model is a scaled-up scenario of low-mass star for- 
, mation: in this case, collimated jets and accretion disk should 
' have the same properties as in low-mass YSOs. The competi- 
tive accretion model suggests that massive stars form in densely 
clustered environments and that disks and collimated jets are 
perturbed by interaction with stellar companions. To distinguish 
between these models, it would be helpful to obtain observa- 
tion of disk/jet systems in a statistical sample of massive YSOs 
to assess the fraction of sources without such system, and to de- 
rive their properties to compare with low-mass YSOs. Unlike the 
still-challenging direct observation of circumstellar disks, out- 
flows/jets extend on larger scales and are easier to observe. 

High-angular resolution studies of low-mass star- forming re- 
gions show that silicon monoxide (SiO) thermal emission is the 
best tool for tracing jets: unlike other tracers such as CO, it 
suffers minimal contamination from swept-up cavities, and it 



allows one to probe obscured regions close to protostars (e.g. 
ICodella et all l2007h . The formation of SiO can be attributed 
to the sputtering of Si atoms from charge d grains in a magne- 
tised C-shock with ve locities >25 km s"' (ISchilke et al.L 119971: 
iGusdorf et al.L 12008 a'). Whether SiO emission originates from 
shocked ambient mater ial or traces the piimary je t itself is still a 
matter of debate (e.g., ICabrit et aU l2007l 1201 2[ and references 
therein). While SiO is extensively observed in low-mass YSOs, 
high-angular resolution studies of SiO thermal emission in high 
mass star-forming regions are still limited to a handful of ob 
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iect s (ICesaroni et al .Ul999^l Huiiter et allll999l:l0iu etani2007t 
Zh ang et all I2007L Zapata e t al.L l2009h with a typical resolution 
of several thousand AU. 

The prominent far-IR source IRAS 17233-3606 is one of the 
best laboratories for studing mas sive star forma tion because of 
its close distance (1 kpc, L eurini et al.L 1201 lah . high luminos- 
ity (~ 1.7 X IO^Lq), and relativel y simple geotnetry. I n a previ- 
ous study in CO with the SMA (ILeurini et all 12009 ', hereafter 
Paper 1), we resolved multiple outflows with high collimation 
factors, extremely high velocity (EHV) emission {\v - ulsrI up 
to 200 km s"'), and flow parameters typical of massive YSOs. 
However, the severe contamination from swept-up components 
and the very elongated beam (~ 5" x 3") did not allow us to 
understand the number of outflows in the region or to identify 
the corresponding driving sources and to assess the occurrence 
of jets. 

In this paper, we present SMA SiO(5^) and complemen- 
tary APEX SiO(5^) and (8-7) observations aimed at separating 
the outflow multiplicity in the cluster We focus on the OFl flow. 
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Table 1. Summary of the observations. 



Centre of the map 

Beam size 
Central frequency 
Velocity resolution 
Bandwidth 



Centre of the map 

Beam size 
Central frequency" 
Velocity resolution" 
Bandwidth 



SMA 

aj2ooo 17''26'"42'-;651 

<5j2ooo - 36°09'38'.'00 

3'.'2 X 2'.'5, P.A. 16° 

217104.980 MHz 

2 km s"' 

2 GHz 



APEX 

aj2ooo 17''26'"42930 

(5j2ooo - 36°09'20'.'00 

3a' 5 

217104.9800 MHz, 347330.5786 MHz 

0.3 km s"', 0.2 km S-' 

IGHz 



Table 2. Summary of relevant positions in the IRAS 17233-3606 
cluster 



mm-continuum 

VLA-2a 

VLA-2b 

VLA-2C 

VLA-2d 

SiO-APEX 

H2-OFI 

SiO-HV" 

SiO-LV'' 

RP 



Q'[J2000] 
Q'[J2000] 
Q'[I2000] 
Q'[J2000] 
Q'[J2000] 
fl'lJ2000] 
Q'[J2000J 
Q'[J2000] 
Q'[J2000] 
fi'[J2000] 



17''26™4Z51 
17'"26"4Z51 
17''26™42^55 
17''26"42':61 
17''26"'42!69 
17''26"4H77 
17''26"'42=:31 
17h26"42':29 
17''26'"4240 
17''26"42U5 



f[J2000] 
'5[J2000] 



f[J2000] 
'5[J20001 



"[J2000] 
'5[J20001 



"[J2000] 
'5[J20001 



"[J2000] 
'^[J20001 



-36°09' 
-36°09' 
-36°09' 
-36°09' 
-36°09' 
-36°09' 
-36°09' 
-36°09' 
-36°09' 
-36°09' 



1 8'.' 32 
18'.' 04 
17'.' 67 
17'.'35 
17'.' 05 
2a.' 00 
14'.'49 
13'.' 74 
15'.' 62 
13'.' 28 



" for SiO(5-4) and (8-7), respectively 



" peak of the HV red-shifted SMA-SiO emission (Fig.lTJ 
'' peak of the LV red-shifted SMA-SiO emission (Fig.[Tll 
'^ peak of the EHV red-shifted CO emission (see Paper I and Fig.O 



which is associated with compact jet-like H2 emission at 2. 12/im 
and extremely high velocity CO, and is possibly driven by a 
zero-age main-sequence (ZAMS) Bl star dZapata et al.Ll2008h . 

2. Observations and data reduction 

2.1. APEX 

We mapped a region of 70" x 70" in the SiO(5^) line at 
217.105 GHz and of 36" x 36" in SiO(8-7) at 347.331 GHz in 
raster mode with the APEX-1 and -2 receivers, respectively, on 
the APE^^Q telescope in May 2010. The spacing between point- 
ing centres was chosen to be 14"at 217 GHz and 9"at 347 GHz. 
The centre of the maps was taken as a[j2ooo] = 17*'26™42';93, 
/6[.T7nnni =-36°0 9'20". The APEX faciUty fast Fourier transform 
spectrometers (iKlein et al.Ll2006l) allow a simultaneous observa- 
tions of 1 GHz bandwidth per receiver, centred on the SiO(5^) 
and (8-7) lines. The velocity resolution of the original data is 
~ 0.3 km s-i at 217 GHz and ~ 0.2 km s"' at 347 GHz. Data 
were converted into Tmb units assuming a forward efficiency 
of 0.95 for both receivers, and a beam efficiency of 0.75 and 
0.73 for the SiO(5^) and (8-7) line, respectively. Pointing was 
checked on the CO(2-l) and (3-2) fines towards RAFGL4211, 
NGC6072 and NGC6302 and was found to have an r.m.s. of 
(3'.'3, 3'.'8) for the APEX-1 data, and of (1'.'4, ^.'6) for the APEX- 
2 data. Data cubes were made with the XY_MAP task in the 
GILDAS softwareQ with a common angular resolution of 30'.'5. 
The final spectral resolution is 3.0 km s ', with an rm.s. of 
0.02 K for SiO(8-7) and 0.1 K for SiO(5^). 

2.2. SMA 

SiO(5-4) fine observations of IRAS 17233-3606 were obtained 
in June 2008 witfi tfie Submillimeter ArrajQ (SMA) in its 
compact-nortfi configuration with baselines ranging in projected 
length from 6 to 74 kA. The central frequency of the lower side- 
band receivers was 217.1049 GHz, while the upper sideband 



' APEX is a collaboration between the Max-Planck-Institut fiir 
Radioastronomie, the European Southern Observatory, and the Onsala 
Space Observatory. 

2 http://www.iram.fr/IRAMFR/GILDAS 

^ The Submillimeter Array (SMA) is a joint project between the 
Smithsonian Astrophysical Observatory and the Academia Sinica 
Institute of Astronomy and Astrophysics, and is funded by the 
Smithsonian Institution and the Academia Sinica. 



central frequency was 227. 1049 GHz. The full bandwidth of the 
SMA correlator was 4 GHz (2 GHz in each band). The SMA 
digital coiTelator was configured to provide a spectral resolution 
of 0.406 MHz (~ 0.56 km s^'). The primary beam at 217 GHz 
has a FWHM diameter of ~ 54". 

The region was covered with a mosaic of two fields 
(arj2000] = 17h26™4Z64, 5[j20()0]=-36°09'18'.'0, ayzooo] = 
17'>26"42'!64, 5[j2()oo]=-36°09'38'.'0). Bandpass calibration was 
made with 3C279. We used Titan for the flux calibration with an 
accuracy of 15-20%, based on the SMA monitoring of quasars. 
The gain calibration was done via frequent observations of 
the quasar 1626-298. The zenith opacity, measured with the 
NRAO tipping radiometer located at the Caltech Submillimeter 
Observatory, was very stable at ~ 0.2, indicating reasonable 
weather conditions during the observations. 

The initial flagging and calibration was performed with the 
IDL superset MIRq The imaging and data analysis was con- 
ducted in MIRIAD (ISault et all Il995 l) and GILDAS. Images 
were produced using natural weighting. The resulting synthe- 
sised beam is 3'.'2 x 2'.'5 with a P.A. of 16°. The final rm.s. cor- 
responds to 0.1 Jy/beam with a spectral resolution of 2 km s"'. 
The conversion factor between Jy/beam and K is 3.4. 

3. Observational results 

3.1. SMA data 

The continuum emission of the region at 217 GHz is shown 
in Fig. [U We performed a 2-D Gaussian fit of the mm con- 
tinuum emission, and determined its peak to be at aj2ooo = 
17h26'"42^.506, Smm = -36°09'18'.'319 with positional eiTors 
on the fit of ~ a.'04. The peak flux is ~ 2.10 + 0.04 Jy beam"', 
the integrated flux 6.69 Jy. The results agr ee very well with th e 
continuum data at 230 GHz presented bv .Leurini et al.l (1201 lal) . 
Figure |2] shows the emission of the SiO(5-4) line ob- 
served with the SMA and integrated over three dififerent veloc- 
ity ranges: the ambient velocity and the blue- and red-shifted 
velocities, which include emission at low (LV) and high veloc- 
ities (HV). The systerri ic velocity of the source is -3.4 km s"' 
(iBronfman et al.L Il996l) . However large uncertainties in the de- 
termination of the ambient velocity were reported in previ- 
ous SMA observations of complex mole cules towards the peak 
of the millimetre continuum emission dLeurini et al.L 1201 Ibl) . 



** The MIR cookbook by Charlie Qi can be found at 
|http://cfa- www.harvard.edu/~cqi/mircook.htmIi 
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Fig. 1. Left: Dashed contours (from lOcr, 0.23 Jy/beam, in step of lOcr) are the 1.4 mm continuum emission observed with the 
SMA. The soUd blue and red contours are the SO(56 - 45) blue- and red-shifted emission from Paper 1 (blue: from 80 to 104 Jy 
beam ' in steps of 10 Jy beam"'; red: 10 and 12 Jy beam '). Right: Integrated emission of the SiO(5^) line: the dashed blue and 
red contours are the LV blue- and red-shifted SiO emission (ubi = [-20, -10] km s"', from 5cr (2.5 Jy/beam) in step of lOcr; v^i - 
[h-10, +39] km s"', from 5cr (4 Jy/beam) in step of lOcr), the thick contour the HV red-shifted emission {v^^ - [+50, +70] km s"'. 
Act (2.4 Jy/beam)). The red dotted contour is the EHV CO(2-l) red-shifted emission (see Fig. 3 of Paper I). The black squares 
mark the peak positions of the EHV CO(2-l) (Rl) red-shifted emission, of the SiO(5-4) LV (SiO-LV) and HV red-shifted emission 
(SiO-HV) , and of the H2 emission in OFl (H2-OFI). In both panels, the grey scale and the solid black contours rep resent the H2 
emission at 2. 12/im. The black crosses mark the positions of the VLA-2a,-2b,-2c and -2d 1.3 cm continuum sources (Zapata et alj, 
I2OO81) : the arrow marks the OFl outflow. The blue and red triangles mark the peaks of the blue- and red-shifted SO(56-45) emission. 



Therefore, we conservatively assumed ambient emission (in 
green in Fig. |2]l from -10 to 10 km s '. Blue- and red-shifted 
emission is clearly detected towards at least two of the out- 
flows identified in the region in Paper 1, OFl and OF2. SiO 
emission is also detected near the H2 emission to the south of 
tiie IRAS 17233-3606 cluster, thus supporting the scenario that 
these H2 structures are also associated with the outflows, as there 
is no H2 emission elsewhere in the map. Ambient emission is 
detected all along OF2 and towards the H2 bow to the north of 
IRAS 17233-3606 . While in our previous study we suggested 
two possible axes for OF2 (Fig.|2] left panel), the SiO(5^) data 
solve the degeneracy and clearly outline its axis (Fig. |2] right 
panel). 

The blue- and red-shifted SiO emission shown in Fig.|2]con- 
centrates along the axis of OFl (see Fig. [TJ. The OFl outflow 
shows interesting properties that deserve a closer look (Fig. [1]). 
The position-velocity diagram of the SiO(5-4) line along the axis 
of OFl is presented in Fig.|3] the SiO emission follows a linear 
velocity-distance relation, typically referred to as the Hubble- 
law, where t he maximum radial velocity is proportional to po- 
sition (e.g., iLada & FichL 11996"). The morphology of OFl as 
traced by SiO is very similar to that pictured by H2. Its red lobe, 
observed in H2 and in EHV CO emission (up to ~ 120 km s '), is 
here confirmed by SiO red-shifted emission at LV and HV veloc- 
ities (up to ~ 55 km s"' at a 3cr level, see Fig.[3]|. Interestingly, 
the HV SiO emission peaks between the apex of the H2 emission 
and the peak of the EHV CO emission. The blue lobe is traced 
by LV and HV CO and SiO (see Fig.[B. The LV SiO emission 
is also very similar to that of SO(56 - 45) (see Fig.[Tl left panel). 



which is the only tracer, together with SiO, that clearly shows 
a bipolar morphology along OFl. Interestingly, the similarity 
between the SiO an d SO emissions follows the recent result of 
Tafalla et aP (1201 Oh . who found that the SiO and SO molecules 
have very similar line profiles in collimated molecular outflows 
driven by Class low-mass sources. Based on the FWHM of 
the LV blue- and red-shifted SiO integrated emission, we infer 
an observed area of ~ 30" squared and a length of ~ 9'.'6, cor- 
responding to a deconvolved area of 22" squared and to a de- 
convolved length of ~ 9'.'0. However, the emission is unresolved 
along the minor axis of OFl . 

It is interesting to compare the SiO and H2 images of the 
OFl outflow with the recent findings towards typ ical low-mass 
protos tars such as HH211 and HH212 (see e.g., i Hirano et al.L 
200^ ICodella et al.L l2007h located at ~ 300-500 pc distance. 
The SiO emission in HH21 1 and HH212 is confined to a highly 
collimated bipolar jet, well detected thanks to sub-arcsec spatial 
resolution of the SMA and Plateau de Bure images. The SiO jet 
can be traced to within 500-1000 AU of the protostar, in a region 
that is heavily obscured in H2 images. On the other hand, the 
brightest H2 knots are located in an outer region and lack a well- 
defined SiO counterpart, probably tracing more powerful shock 
where SiO does not form/survive or where the SiO excitation is 
extremely high and consequently the low-J mm-transitions are 
weak. In addition, the position-velocity diagrams suggest that 
the SiO emission increases in velocity with increasing distance 
from the driving source. The present OFl images, although ob- 
served with a poorer spatial resolution, seem to suggest a similar 
picture: (i) the position of the SiO red lobe is not correlated with 
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Fig. 2. Left panel: Integrated emission of the CO(2-l) blue- and red-shifted wings. The solid contours show the EHV blue- and red- 
shifted emission; the dashed contours mark the HV blue- and red-shifted emission (see Paper I, Fig. 3). The dashed and solid lines 
outline the possible molecular outflows (OFl, OF2, OF3); the dotted line marks the alternative direction of the 0F2 outflow. The 
central box is the region mapped in Fig. [T] Middle panel: Integrated emission of the SiO(5-4) line: the dashed blue and red contours 
are the LV blue- and red-shifted emission (ubi = [-20, -10]km s"', from5cr(2.5 Jy/beam)in steps of lOcr; iifd - [-1-10, -1-39] km s"\ 
from 5cr (6 Jy/beam) in steps of lOcr), the thick contours the HV blue- and red-shifted emission (^bi - [-50, -20] km s"', from 5cr 
(4 Jy/beam) in steps of lOcr; v^^ - [-(-50, +70] km s"', Act (2.4 Jy/beam)). The green contours are the integrated emission at ambient 
velocities (iig, = [-10, +10] km s"' ; from 5cr (3 Jy/beam) in steps of lOcr). Right panel: Integrated emission of the SiO(5-4) line in 
the velocity ranges used for modelling the SiO emission in Sect. ID the blue and red contours are the blue- and red-shifted emission 
in the ranges (ubi - [-30, -20] km s"', from 5cr (4 Jy/beam) in steps of lOcr; v^^ - [+10, +39] km s"', from 5cr (6.5 Jy/beam) in 
steps of lOcr). In the three panels, the solid line marks the axis of OF2 as identified in SiO. The arrowed line marks that of OFl . The 
solid circle represents the APEX 30'.'5 beam centred on a[y2000] - 17^26'M1!77, ^[72000] = -36°09'20", the position used in Sect.|4] 
for our modelling of the SiO emission (SiO- APEX, see Table|2|i and here labelled by a black triangle. In both panels, the grey scale 
represents the H2 emission at 2.12jum (Paper I). 



the brightest H2 knot, and (ii) the SiO-LV emission peaks closer 
to the driving source than the SiO-HV peak. Only future higher 
spatial resolution observations, possibly with ALMA, will allow 
a proper characterisation of the region. 

Comparisons between the SMA SiO(5-4) flux and the flux 
observed with the APEX telescope in the same transition are 
discussed in the next section. 



3.2. APEX data 

Figure |4] shows the SiO(5^) and (8-7) maps from APEX at 
their original angular resolution, 30'.'5 forSiO(5-4) and 19'.' 1 for 
SiO(8-7). Given the low resolution of the APEX data, the emis- 
sion from the different outflows in the region is not resolved. 
When averaged to the common resolution of 30'.'5, both SiO(8- 
7) and (5-4) maps show a maximum at Q'[y20()0] = 17^26™4P!77, 
'5[y2()00] - -36°09'20", the centre of the circle (representing the 
APEX beam) partially seen in Fig.[2](right panel). This position 
has no counterpart in other tracers: it is situated at ~ (-9", -2") 
from VLA-2a, the strongest cm source of the region, and at 
~ (-8", -4") from the centre of the LV SiO red-shifted emission 
in OFl. We believe, however, that the offset between SiO SMA 
and APEX data is due to the poor sampling of the APEX map. 



In Fig.|5]we present the SiO(5^) emission observed with APEX 
in three different velocity ranges (total emission, red-shifted and 
blue-shifted emission). The map clearly shows that SiO is spread 
over a large area and is associated with several shocks from the 
IRAS 17233-3606 cluster detected in H2. 

The APEX spectra extracted at the peak of the single-dish 
SiO emission are shown in Fig. |6] Although the APEX beam 
at this position covers emission close to the ambient velocity 
from OFl and OF2, the red- and the blue-shifted emission is 
dominated by OFl both spatially and in intensity, as discussed in 
Sect. l3.3l Therefore, in the following we assume that the spectra 
shown in Fig. |6] are associated with OFl only and represent its 
integrated spectra over the 30'.'5 APEX beam. 

The SiO(8-7) line shows bright wings extending to 
+40 km s"', i.e. velocities slightly higher than those of the 
SiO(5^) transition. Fig.|6]also shows that the SiO(8-7)/SiO(5- 
4) intensity ratio remains almost constant (~ 0.8) at ambient ve- 
locities while it increases to values around 1.7 at red-shifted ve- 
locities, and 2.2 in the blue wing. If we assume that both SiO 
lines are tracing the same gas, this result suggests an increase 
of excitation conditions (discussed in Sect. |4| at high velocities, 
similar to the results obtained towards the two low-mass Class 
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Fig. 3. Position-velocity diagram of the SiO(5^) transition 
(SMA data) along the OFl axis shown in Fig. [T] The vertical 
line marks the mm continuum p eak, the ho rizontal line the am- 
bient velocity (-3.4 km s~'. Bronfman et al.i 1996) . Contours are 
from 2cr (0.2 Jy/beam) in steps of 4cr. 



sources L1448-mm and LI 157-mm (iNisini et al.L 120071) through 
SiO observations. 



3.3. Comparison between the SMA and the APEX SiO 
emission 

We compared the SiO(5^) APEX spectrum at the centre of 
the single-dish map with the corresponding SMA spectrum inte- 
grated over the APEX beam. The SMA total flux corresponds to 
41% of the single-dish emission, 57% in the blue-wings (which 
dominate the spectrum at this position), and ~ 34% at ambi- 
ent velocities. Given the broad band of the APEX and of the 
SMA receivers and the richness of the molecular spectrum of 
IRAS 17233-3606 dLeurini et al.L llool and Paper I), several 
spectral features are detected simultaneously to SiO(5^) in both 
datasets. For emission in the band from complex molecules, 
most likely probing the hot core, which is expected to be com- 
pact, fluxes between 80 and 100% are recovered, hence confirm- 
ing that the difference between the SMA and APEX SiO spectra 
is due to missing flux and not to calibration issues. 

We also compared the SMA spectrum integrated over the 
APEX beam and centred on the peak of the APEX SiO emis- 
sion with the corresponding APEX spectrum. The overlay of 
the two spectra is shown in Fig. Q The SMA recovers 98% 
and 68% of the APEX flux in the red ([-h10,-h39] km s"') and 
blue low-velocity ([-20,-10] km s"') ranges plotted in Fig. |2j5. 
In the APEX spectrum, emission at high-velocity range (t^id - 
[h-50, -h70] km s"' and vu < -30 km s"') is not detected. 
However, emission is detected at high velocities in the [-30,- 
20] km s"' range. In this range of velocities, the SMA recovers 
the whole APEX flux (Fig. |7]). Therefore, in the range of veloc- 
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Fig. 4. SiO(5^) (top panel) and (8-7) (bottom panel) maps ob- 
served with the APEX telescope. The velocity axis extendes 
from -95 to 85 km s"' for both data sets, while the temperature 
scale is covers -0.4 to 3 K. The data are shown with the origi- 
nal angular resolution of 30'.'5 (SiO(5-4)) and 19'.' 1 (SiO(8-7)). 
The centre of the maps is given in Table [1] The black triangle 
marks the peak position (SiO- APEX, see Table |2]i of the inte- 
grated intensity SiO(5^) and (8-7) maps (in the velocity range 
-40 < y < 40 km s ') smoothed to the angular resolution of the 
SiO(5-4) data, 30'.'5. 



ities [-HlO,+39] km s"' and [-30,-20] km s"', the SMA data are 
not afifected by filtering of large structures and the morphology 
of Fig.|2}; (which shows the SMA SiO integrated emission in the 
ranges [-i-10,+39] km s"' and [-30,-20] km s"') is reliable. We 
show in Appendix |A] a joint deconvolution of the inteferome- 
ter and single-dish data performed on integrated intensity maps. 
However, combining the inteferometer data with the single-dish 
map is not needed for our purposes because the modelling of 
the SiO emission (see Sect. |4]l is limited to the velocity ranges 
[-1-10,-1-39] km s"' and [-30,-20] km s"' where the SMA recovers 
the whole flux observed with APEX. 

To estimate the contribution of OFl to the APEX spectrum at 
the peak of the SiO emission, we computed the SMA fluxes inte- 
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Fig. 5. Total integrated intensity (-40 < y < 40 km s"'; grey 
shades are from 10% of the peak emission, 43 K k ms"', in 
steps of 5%) of the SiO(5^) transition observed with APEX. 
The red and blue contours are the red- (5 < y < 40 km s ') 
and blue-shifted (-40 < u < -5 km s"') integrated emission 
observed with APEX (contours are 50% of the peak emission 
(9.4 K k ms"' for the red emission, 17.6 K k ms"' for the blue 
emission) in steps of 10%). The green contours represent the 
2.12jum H2 emission. The white arrows label the OFl and OF2 
outflows. The black circle marks the APEX beam centred on the 
peak of the APEX SiO emission, as in Fig.|2] 



grated over the APEX beam centred on (^[yaooo] - 17''26™41'577, 
<5[y2000] = -36°09'20" and over a region enclosing OFl. The 
blue-shifted ([-30, -20] km s"' ) SMA flux associated with OFl 
is 66 Jy, corresponding to 66% of the total blue-shifted flux over 
the APEX beam. Similarly, the red-shifted ([-hlO, -i-39] km s^^) 
SMA flux in OFl is 71Jy (53% of the total flux over the APEX 
beam). Therefore, OFl contributes up to ~ 60% to the integrated 
SiO emission observed with APEX in the velocity ranges used 
in Sect.|4]for the modeUing of the SiO spectra. 



4. SiO modelling 

In principle, specific modelling complications arise in the con- 
text of high-mass star formation, where the strong radiation field 
of the central object may affect both the physical and chemical 
structure of the shock and the excitation of molecules in the out- 
flow region (such effects on the shock models are currently being 
studied, but their influence has not been properly assessed yet). 
The potential multiplicity of the sources or their outflows is also 
likely to generate modelling difficulties. Owing to the lack of in- 
formation on the region, and as a first step towards a more com- 
plete description, we decided to model our SiO observations in 
two ways that are classically used to describe the SiO emission 
arising from shock re gions generated in outflows of low-mas s 
forming stai's (e.g., Nisin i et al.l 120071: iGusdorfetall 1201 li) . 
The first one, hereafter referred to as 'shock-LVG', is described 
in Sect. 14.11 and is based on the combination of a shock code 
and a radiative transfer code based on the large-velocity-gradient 
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Fig. 6. Upper panel: APEX spectra of the SiO(5^) and (8-7) 
lines in the OFl region around IRAS 17233-3606 observed at 
ffpooo] = 17h26"4P:77, dumm] = -36°09'20". Lower panel: 
SiO(8-7)/SiO(5^) line ratio as function of velocity (black line). 
The red line is the SiO(8-7) spectrum. The dashed lines mark 
the velocity ranges used in Sect. |4] blue [-30, -20], and red 
[+10,+39]kms-i. 
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Fig. 7. SiO(5^) spectra from APEX (black fine) and the SMA 
(red line) observed towards OFl at ajyiooo] = 17''26'"4P.77, 
^[72000] = -36°09'20". The dashed lines mark the velocity 
ranges used in Sect. |4] for modelling the SiO emission, blue 
[-30, -20], and red [-hIO, -h39] km s"'. 



(LVG) approximation. The second one (Sect. I4.2l i. hereafter re- 
ferred to as 'slab-LVG', consists of a simpler LVG description. 
In the first one, the shock is approximated by a layer of points 
with varying physical and chemical conditions from one point 
to another In the second, the usual LVG physical parameters are 
given for one point that epitomises the whole layer of the shock- 
LVG method. Although still distant from the probably complex 
'real' configuration of multiple shocks at work within the APEX 
beam, we believe that this two-sided approach, especially the 
shock-modelling one, is the most rigorous so far available to 
constrain the abundance of SiO and link it to the mass of the 
associated outflow. Indeed, from a computational point of view, 
it is still impossible to design a model that would account for the 
precise geometrical, dynamical, physical, and chemical structure 
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of the region. Moreover, no observational facility can currently 
provide the very high-angular resolution necessary to separate, 
for instance, a 'wind shock' from an 'ambient shock' component 
identifi ed in the idealistic models presented, e.g., bv iHollenbachI 
(IT997h . 



4.1. Shock-LVG analysis 

We first base our modelling of the SiO emissi on on the C- 
type shocks interpretation presented in details in iGusdo rf et al.l 
(2008a). The SiO emission is accounted for by the combination 
of two codes: a state-of-the-art shock model and a radiative- 
transfer code based on the LVG approximation. The t w o code s 
codes are presented in detail in iGusdorf et al.l (l2008alfbl 1201 ih . 
We nevertheless describe the main elements relevant to our study 
in the following paragraph. 

The shock code calculates the dynamical, physical, and 
chemical structure of one-dimensional, stationary C-type shock 
models. The fractional abundances of over 125 species, linked 
by over 1000 reactions, are computed at each point of the shock. 
The gas-phase abundance of H2 is calculated, base d mostly 
on the treatment described in iLe Bourlot et al.l (l2002h . The lo- 
cal abundance of H2 is of particular importance in the context 
of SiO modelling, as it is the most abundant species to colli- 
sionally excite SiO. The excitation of its own 150 first levels 
(the most likely to be populated in the considered shocks) is 
also provide d, making use of the H -H2 collisional rate coef- 
ficients from 'Wrathmall et al.' ('2007*). Specific grain processes, 
listed in Flower & Pineau des Forets (2003), are taken into ac- 
count. In particular, in C-shock models, the grains undergo the 
sputtering impact of neutral particles at the speed of the am- 
bipolar diffusion. Because the grain cores are partly made from 
silicon-bearing material, this sputtering leads to the release of 
elemental Si in th e gas phase (Caselli et al. 1997; Schilke et al. 
[19971; iMav et al.1 (2000), where it reacts with O2 and OH to 
form SiO. In the shock code, the rates of these two reactions 
are identical and correspond to the low-tempe rature measure- 
ments of the Si + O2 — > SiO -1- O reaction by iLe Picard et al.l 
(I2OOII) . Less prominent routes to the formation of SiO, such as 
Si+(OH,H)SiO+(H2,H)HSiO+(e,H)SiO, are also included in our 
shock code. 

The magnetic field parameter, b, defined by 

B(//G) - b X ^J^^l(cn^^ where hh is the pre-shock 
density, is set to 1 . The shock velocity ranges from the minimum 
value likely to generate sputtering of the grain cores (over 
25 km s"'), to the maximum value permitted by the peculiar 
physics of those sh ocks and by the choice of such parameters 
dFlower et al.l l2003l) . The age of the shock region can also be 
varied in our models, by varying their final flow time, that is the 
necessary time for a particle to flow from the pre-shock to the 
post-shock region. In practice, the final flow time is the time at 
which we stop the calculations. 

The necessary outputs of the shock model, such as the tem- 
perature profile and fractional abundances of relevant species 
(SiO and colliding species) are then used as inputs of the 
radiative-transfer code, which also makes use of SiO spectro- 
scopic para meters (Einstein co efficients, collisional rate coeffi - 
cients) from lDavou & Balanqal ([2006) and lSchoier et al.1 (l2005l) . 
and provided by the LAMDA databasqfl to solve the statistical 
equilibrium equations and calculate the level populations at each 
point of the shock. The set of equations is solved thanks to the 
escape probability formahsm. The fine emissivities and optical 
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depth are also computed at each point of the shock layer for all 
considered transitions. Finally, we can compare the integrated 
intensity over the whole shock region, as resulting from the mod- 
els, with the observations. 

The 30'.'5 beam of the APEX observations encompasses the 
blue- and red-lobe of OFl (see Figs. |2] and |5] and also Fig. |6] 
for our definition of the velocity intervals associated with these 
lobes), which correspond to different shock regions detected at 
high velocities. It also comprises an 'ambient velocity compo- 
nent' that we interpret as a consequence of the propagation of 
the shock wave. Indeed, this narrow, ambient-velocity compo- 
nent has already been seen in bipolar outflows associated to 
low-mass forming stars, and it could arise either from deceler- 
ating Si-enriched post-shock gas, mixed with the ambient gas 
(e.g., iLefloch et al.lll998tlCodella et all 1 19991) . or from a slow 
reverse C-type shock brough t to rest by the much denser am- 
bient medium (iGusdorf et al.L l2008al) . In any case, this part of 
the SiO emission is not covered by our shock modelling, which 
solely describes the more dynamical SiO emission detected in 
the blue- and red-wings of our observations. Owing to its asso- 
ciated velocity, this SiO emission is interpreted as arising from a 
shock of outflow- or jet-like origin, a view that is also supported 
by the maps shown in Figs.[T||2]and|5] 

We hence compare the red-lobe and blue-lobe integrated 
intensities of the SiO(5-4) and (8-7) lines with the results 
generated by the grid of C-sh ock models previously intro- 
duced in IGusdorf et al.l (l2008al) . for pre-shock density values 
of lO'*, 10^, and 10^ cm"^. The observed integrated intensities 
are 2.1 + 0.4 and 4.4 + 0.9 K km s"' in the blue wings and 
3.8 + 0.8 and 6.7 + 1.3 K km s^' in the red wings, for SiO(5- 
4) and (8-7) respectively. The velocity ranges used in the anal- 
ysis are [-30, -20] km s ' for the blue-shifted emission, and 
[-1-10,-1-39] km s ' for the red-shifted emission as defined in 
Sect. l3.3l and labelled in Fig.|6] The error bars associated to these 
values correspond to the 20% calibration uncertainty intrinsic to 
our measurements. 

From our SMA maps, we derived an observed area of ~ 
30 arcsec^ for the SiO(5-4) emission in OFl, which corresponds 
to a deconvolved area of ~ 22 arcsec-, although the emission 
is unresolved along the minor axis of the emission. Therefore, 
we adopt this value as an upper limit to the size of the SiO 
emission associated with OFl. A second estimate of the size 
of OFl comes from the H2 emission. The assumption that SiO 
and H2 have similar sizes in protostellar outflows is justified by 
high-angular resolu t ion observations of lo w-mass systems (e.g., 
iHirano et all l2006t ICodella et all l2007l) . which clearly show 
that SiO knots have near-infrared H2 counterparts except at high 
extinctions associated with dense protostellar envelope where 
H2 is not detected. For the H2 counterpart of the red lobe of 
the SiO(5^) emission shown in Fig.[T] we measure an emitting 
size of ~ 6 arcsec-. For the blue lobe, no H2 emission is de- 
tected, possibly owing to extinction. However, because OFl is 
seen symmetric in both SiO and SO (Fig.lTJ, we assume an emit- 
ting size of ~ 6 arcsec^ for the blue-lobe, resulting in a total area 
~ 12 arcsec^. Finally, we assume that the size of the emitting 
region is the same for both SiO lines; hence, the beam-filling 
factor is the same for both transitions, since the SiO(8-7) data 
were smoothed to the resolution of the SiO(5^) observations 
(Sect. lZTt . This minimises possible biases due to beam dilution. 

We base our selection of a 'satisfyingly fitting model' on 
two criteria. First, we use the integrated intensity ratio of the 
SiO(8-7) to th e SiO (5^) lines. Ac cording to Fig. [8] (adapted 
from Fig. 9 of IGusdorf et al.ll2008al) . the only model reproduc- 
ing the observations is that with riu = 10^ cm"^, b - I, and 
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t^s = 32 km s ' . Below this velocity, the sputtering and excitation 
conditions are not efficient enough, whereas above it, too much 
SiO(8-7) emission is generated to match the observations. Then, 
we check that the maximum line temperatures predicted by the 
model agree with the observed ones, which can be read in main 
beam units in Fig. |6] around 0.6 K in the red and blue wings, 
with a maximum slightly higher for SiO (8-7) than for SiO (5- 
4). Assuming that ~60% of the APEX SiO emission originates 
in OFl (Sect. [331 l. one must then account for a main beam max- 
imum temperature of ~0.36 K. To convert this into brightness 
temperatures unit, we use both the upper limit of 22 arcsec^ and 
the H2 emission size of 6 arcsec^ for each lobe of SiO emission. 
The beam size of the APEX spectra being BCK'S for both transi- 
tions, the 0.36 K main beam value amounts to ~ 12 and ~44 K af- 
ter correcting for dilution. This latter value agrees very decently 
with our models predictions (Fig.|8]l, where the line temperature 
profiles are displayed over the whole layer of the shock for the 
(5^) and (8-7) transitions, showing maximum temperatures of 
respectively 49 K and 51 K. We note that such a satisfactory 
agreement cannot be found for other models of our grid. This 
seems also to suggest an emitting size for the SiO(5^) compa- 
rable to that of H2, ~ 6 arcsec^ per lobe. 

We can then extract the surface density corresponding to this 
best model, and derive two values corresponding to 500 and 
1000 years, which are the typical upper limits for the age of our 
shocks given in Paper I. We obtain 0.12 and 0.25 g cm"^. The 
second value is higher than the first one because an older shock 
age implies a more-extended post-shock region, where the gas 
is denser. We first derive the total outflow-mass comprised in 
the upper limit of our SiO emitting region size, ~ 22 arcsec^, 
adopting a distance of 1 kpc. We obtain a total mass of 0.3 Mq 
with the 500 years limit, and 0.6 M© with the 1000 years limit. 
If we alternatively use the total emission region size of H2, 
namely ~ 12 arcsec^, we obtain a corresponding total outflow- 
mass of 0.2 M0 for 500 yr, 0.3 M© for 1000 yr approximately. 
The calculated masses agree with those inferred in Paper I (0.6- 
0.9 Mq) and imply a luminosity of the driving source of OFl 
of ~ lO-' Lq based on the empirical outflow mass-luminosity 
relation of Lopez-Sepulcre et al. (2009). This luminosity corre- 
sponds to that of a B3-type ZAMS star (e.g., lPanagialll973h . 
However, our estimate of the total mass of the outflow has to 
be considered as a lower limit since only a limited range of ve- 
locities was used in the analysis. The inferred luminosity of the 
powering source of OFl is also probably a lower limit because 
the o utflow mass-luminosity relation of iLopez-Sepulcre et alJ 
(I2009.) is based on CO measurements. 



4.2. Slab-LVG analysis 

To verify these results and to constrain the excitation conditions 
of the OFl SiO outflow, we also performed a 'slab' analysis of 
the OFl region using only an LVG code. We investigated the 
column density of SiO (A^sio) in the range 10" - 10'^ cm"^, H2 
densities (hhj) of lO-' - 10*^ cm""* and kinetic temperatures (Tkin) 
of 10-600 K. Figure [To] shows our typical model results for a ki- 
netic temperature of 100 and 500 K. We analysed the line inten- 
sities at LSR velocity as well as at the average velocities of the 
blue- (-25 km s"') and red-shifted (-1-15 km s"') lobes. Figure 
[TO]indicates temperatures higher than 100 K for the red and blue 
lobes: lower temperatures would imply unreasonable densities 
higher than 10^ cm"^. The tin^ values for r^n = 500 K are ^ 5 
10^ cm"^. Since this is very likely an upper limit to the kinetic 
temperature of the gas, these densities have to be considered as 
lower limits. The inferred SiO column densities are ~ 0.5-1 10'"* 
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Fig. 8. Integrated intensities of the SiO lines relative to SiO(5- 
4). The horizontal lines show the values of the SiO(8-7) to (5-4) 
ratio, whether it is the one in the blue or red lobe of the outflow. 
The vertical line marks the position of the (8-7) transition on 
the X-axis. The colour curves provide the results for the fraction 
of our C-type shock models grid relevant to our study. For all 
these models, the pre-shock density is nu = 10'' cm"^, and the 
magnetic field parameter is /? = 1. Only the shock velocity u^ 
varies, and is indicated for each model on the corresponding line. 
The blue- and red-velocity range are defined in Fig. |6] 



cm for both lobes. On the other hand, at ambient velocity we 
can infer A^sio of 2 - 5 10'^ cm"^ and densities higher than 10^ 
cm"^ without any clear constraint on kinetic temperature. 

For consistency purposes, we cross-checked the results 
yielded by our two modelling methods, shock-LVG vs. slab- 
LVG. To do so, one must return to the shock models and extract 
the ranges of variation over the shock layer of the parameters, 
which can only then be compared with the single-point values 
of the same parameters, which prevail in the slab-LVG analysis. 
The shock code outputs consist of a series of profiles of phys- 
ical and chemical quantities (e.g., temperature, abundance) cal- 
culated over the shock layer Our LVG code, when used in com- 
bination with the shock code, generates a series of outputs in 
the form of similar profiles of emissivities quantities (e.g., level 
populations, line emissivities). For the model that was chosen to 
describe the OFl observations, Fig.|9]shows the neutral tempera- 
ture and the local emissivity of the SiO (5^) line over the whole 
shock layer, which is represented on the x-axis by means of the 
neutral flow time of the shock, i.e. the necessary time for a neu- 
tral particle to fly through the shock region. This representation 
allows us to empirically define the SiO (5^) emission region 
as the region where the emissivity of the line is significant. By 
significant we mean that the emission region consists of all the 
points where the emissivity of the line exceeds 1 % of the max- 
imum line temperature, that is about 0.5 K for the considered 
model. We then compare the range of variations over this emis- 
sion region, of the quantities that are also used as input parame- 
ters in the slab-LVG method. This way, we performed a proper 
comparison between the range-values inferred from the shock- 
LVG analysis and the single-values yielded by the slab-LVG 
method. These are reported in Table |3] In our shock-LVG mod- 
els, only one parameter is not expressed in the form of a range 
of values: the SiO column density because it corresponds to the 
SiO density integrated over the whole shock layer Note that the 
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Fig. 9. Neutral temperature (black line) and the SiO (5-4) and 
(8-7) line temperature (red and blue lines) profiles over the 
whole shock layer, seen face-on, versus the neutral flow time 
for our best-fitting model. The characteristics of this model are 
indicated in the figure, as well as the location of pre- and post- 
shock regions. Our definition of the SiO (5^) emission region 
(Sect. 14.21 ) is also delimited by the two red open circles and is 
represented by the horizontal dashed red line. In this diagram, 
the line of sight lies on the plane of the figure, in the direction 
that goes from the pre- to the post-shock region. 
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Fig. 10. Volume density ne, versus SiO column density, A^(SiO) 
for slab-LVG models for fii„ = 100 K (left panel) and 500 
K (right). Each continuous curve corresponds to the observed 
SiO(8-7)/SiO(5^) intensity ratio (after convolution to the same 
resolution, see Sect. 12.1b . as measured at ambient velocity (- 
3.4 km s"', green curve) and at typical red- (-1-25.0 km s"', red 
curve) and blue-shifted (-15.0 km s ', blue curve) velocities. 
Dashed lines are for the ratio uncertainties. The dotted curves 
correspond to the APEX SiO(8-7) values at the same velocities, 
after correction for beam dilution using the OFl size as derived 
from the SMA data. The corresponding uncertainties are negli- 
gible and are not shown for the sake of clarity. The symbols (cir- 
cle: ambient; square: red; triangles: blue) indicate the allowed 
solutions. 



definition of the LVG parameter is the same in both the slab- 
LVG and shock-LVG cases, but that its expression differs. For 
the slab-LVG modelling, it is simply equal to A^sio (cm"^) / Av 
(km s"'), which is equivalent to the shock context formulation 



w 10 X «(SiO)(cm ) / grad_!; (s ), where grad_i; is the veloc- 



ity gradient, and where the 10^ factor only arises because of the 
conversion from km to cm. 

The results shown in Table [3] are remarkably consistent and 
seem to support the view that tools that have been proved effi- 
cient in the context of low-mass star formation can be used in 
higher mass cases, with slightly increased densities. However, 
our best-fitting shock model was only approaching the observed 
SiO line ratios, and despite additional investigations in the pa- 
rameter space, we were unable find a model likely to match 
(8-7) / (5^) ratio values as high as 2, for instance. If con- 
firmed as characteristic of massive star formation processes, 
such a high ratio could be accounted for by considering grain- 
grain interactions within the shock layer Such interactions have 
been described as important at high pre-shock densities, and 
their feedback effects on both the chemical and dynamical struc- 
ture of shocks have been studied by Guillet et al. (2007, 2003, 
'2011). Unfortunately, those authors could not generate any self- 
consistent solution for the pre-shock densities at work here, 
namely for «h - 10^ cm"^. Alternatively, the influence of a 
strong radiation field generated by the driving source of the out- 
flow both on the shock structure itself and on the correspond- 
ing molecular excitation might yield high SiO line ratios. The 
rather complex implementation of such processes in shock mod- 
els is work in progress. The present analysis indicates volume 
densities higher than lO*" cm""* (see Table 3). lNisini et al.l (120071) 
performed an LVG analysis of SiO emission towards the nearby 
prototypical Class low-mass objects L1148 and L1157, using 
single-dish data on angular scales of 10". No tight constraints 
have been obtained for kinetic temperature, whereas the vol- 
ume densities we re found to be between 10^ and 10^ cm"-'. 
In addition, Gusdorf et al.l (12008 a) modelled the SiO emission 
from Ll 157 using a C-shock code leading to pre-shocked densi- 
ties of lO'^-IO^ cm"-'. Thus, although completely different spa- 
tial scales are involved, the volume densities here inferred for 
IRAS 17233-3606 appear to be consistent with those derived for 
L1448andL1157. 

Although the results from the slab LVG approach are consis- 
tent with those inferred through the more rigorous shock-LVG 
analysis, we stress here the importance of shock models for de- 
riving physical parameters such as the SiO abundance and there- 
fore the mass of the outflow: while in the slab-LVG method one 
has to rely on assumed values of Xsio, ^sio is consistently com- 
puted in the shock model, thus decreasing the uncertainties in 
the estimate of the outflow mass. 



5. The nature of OFl 

The present data sets trace the outflow activity driven by the 
YSOs in the cluster, and contribute to draw a clearer picture of 
the region than that obtained by our previous CO data (Paper I). 
Our previous observations of OFl were limited to EHV red- 
shifted CO(2-l) emission and to the detection of H2 emission 
also constrained to only the red lobe of OFl . The bipolar nature 
of OFl was suggested only by the SO data and by water maser 
spots. Our present SMA data confirm that the emission associ- 
ated with OFl is due to a bipolar collimated outflow well traced 
by SiO in both red- and blue-shifted emission. This is consistent 
with findings for jets driven by Sun-like protostars, where it is 
possible to see that CO traces not only the jet but also the walls 
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Table 3. Properties of SiO-emission regions inferred from face-on C-type shock models used in combination with an LVG code, 
and from homogeneous slab-LVG models. 



medium 
properties 



red lobe 



slab-LVG 

blue lobe 



ambient 



shock-LVG 

best model 



Tkin (K) 


100-500 


100-500 


n.c. 


15-1745 


nn, (cm"') 


> 5 X 10"^ 


> 5 X 10« 


> 10^ 


3.3x10^-2.4x10' 


A'sio (cm"-) 


5x10'^^ 


10" 


5x10" 


2.4X10''* 


LVG parameter 


2.5xl0'2 


5xl0'2 


2.5xl0i2 


6.4xl0"-5.0xl0'^ 



of the cavit y opened by the je t itself (e.g.- lLee et al.ll2007l their 
Fig. 5, and ICabrit etal]|20 1 2l their Fig. 1). On the other hand, 
SiO emission is usually associated with primary jets in the case 
of low-mass YSOs (e.g., Hir ano et al., 2006; Cabrit et al., 201.1) 
and generally suffers minimal contamination from infalling en- 
velopes or swept-up cavities. 

Our data strongly suggest that the emission associated with 
OFl traces a bipolar jet from the IRAS 17233-3606 cluster This 
view is additionally supported by the high velocities detected in 
SiO, and by the collimated morphology of OFl, which is ~ 6 
according to the H2 images, among the larges t ever reported in 
high-mass star formation (see for exam ple Beuther et al.l l2002i 
l2004tlSu et ani2004tlWang et al.Ll201 Ih . This value corresponds 
to the ratio between the major semi-axis and the minor semi- 
axis of OFl as derived from the H2 map. Additionally, the high 
temperature (T > 100 K) derived through the SiO analysis is 
suggestive of material very close to the primary jet. 

Although the driving source of OFl cannot be identified 
yet given the still relatively low resolution of the SMA data 
(~ 3000 AU), the best candidate source is VLA-2a, which is the 
strongest cm-continuum source in the region and also the com- 
pact source closest to the H2 jet-like feature. Alternatively, VLA- 
2b could also be a possible candidate for the powering source of 
OFl . In principle, we cannot rule out the hypothesis that the driv- 
ing source of the observed SiO bipolar emission is a low-mass 
YSO in the cluster. However, there is some indirect evidence that 
OFl is indeed driven by a massive YSO. First, we determined the 
total mass of OFl through the shock modelling of SiO emission 
without any assumption of the SiO abundance (Sect. l4Tb and us- 
ing the area inferred from the H2 emission. From this mass, we 
inferred a luminosity of > 10^ Lq for the OFl driving source, 
which agrees with an earlier than B3-type star The continuum 
flux of VLA-2a in the cm (Zapata et al.. 2008) is also compati- 
ble with t hat of an early B ZA MS star. The cm continuum flux of 
VLA-2b (IZapata et al.L 120081) is also compatible with an early B- 
type star, although less massive than VLA-2a (a B3 ZAMS star 
compared to a Bl ZAMS star). To our knowledge, this is one of 
the first indirect evidences of a SiO collimated bipolar outflow 
from an early B type YSO, although only higher spatial reso- 
lution maps will allow a direct proof of its existence. Previous 
studies reported red- and b lue-shifted SiO emiss i on from mas- 
sive star- formin g regions (jC esaroni et al., 1999; Hunter et al., 
[T99lloiu"erani2007tlZhang et al.. .2007: .Codella et al .2013.) 
Among these, IRAS20126 -H4104 (e.g., ICesaroni et all Il999t 
ICaratti o Garatti et al.L 120081) probably represents the best exam- 
ple of a jet from a massive YSO, since the nature of the jet is 
proved by SiO, H2 and [Fe 11] obs ervations. However, t he mass 
of the powering source is ~ 7 Mm ([Cesaroni et al.Ll2005h . In one 
case (IRAS 18264-1 152. l0iu et all l2007h SiO is also detected 
at high velocities, as in the present case of OFl. However, the 
difference between these studies and our case is that several fea- 
tures typical of low-mass jets (collimation, emission at high and 
extremely-high velocity, the parameters of the gas derived by 



the shock model) are all found in OFl while the other examples 
show only some of these properties. 

6. Conclusions 

We reported and analysed SMA and APEX SiO observations of 
the IRAS 17233-3606 star-forming region. The SMA data al- 
lowed us to clearly resolve the OFl outflow with a resolution of 
~ 3000 AU and identify its SiO emission as due to a collimated 
bipolar outflow associated with EHV CO, H2 and SO. 

The OFl SiO(5-4) and (8-7) APEX specti-a suggest an in- 
crease of the excitation conditions at high velocities and point to 
hot gas close to the primary jet, as found in outflows from low- 
mass YSOs. From the shock-LVG analysis of the APEX spec- 
tra, we derived a mass of > 0.3 Mq for OFl, which implies a 
luminosity L> 10^ L© for the OFl driving source. Although the 
driving source of OFl cannot be yet identified at the resolution 
of our current SMA data, the best candidate source, VLA-2a, has 
a continuum flux in the cm that also agrees with that of a ZAMS 
star with a luminosity of a few lO-' L©. 

The high angular resolution spectroscopic observations that 
we collected on OFl in IRAS 17233-3606 support a picture sim- 
ilar to that seen in studies of jets and outflows from low-mass 
YSOs. The SiO emission at LV peaks closer to the driving source 
than the SiO-HV peak; moreover, SiO traces a region of the out- 
flow/jet system heavily obscured in the H2 image. Finally, the 
SiO and SO molecules have a very similar distribution as found 
in collimated outflows driven by low-mass Class YSOs. 

In terms of shock modelling, we found it possible to model 
SiO emission in shock regions generated by high-mass proto- 
stars with the same tools developed for shocks from low-mass 
YSOs. The shock parameters inferred by the SiO modelling are 
compatible with those found for low-mass protostars, only with 
higher pre-shock densities. The same conclusions arise from a 
simpler LVG slab analysis. 

Acknowledgements. The authors would like to thank the anonymous referee 
who helped to clarify the results and improve the discussion. A. Gusdoif ac- 
knowledges support by the grant ANR-09-BLAN-0231-01 from the French 
Agence Nationale de la Recherche as part of the SCHISM project. 



References 

Beuther, H., Schilke, P., & Gueth, F. 2004, ApJ, 608, 330 

Beuther, H., Schilke, P., Sridharan, T. K., et al. 2002, A&A, 383, 892 

Bonnell, I. A., Lai'son, R. B., & Zinnecker, H. 2007, Protostars and Planets V, 

149 
Bronfman, L., Nyman, L.-A., & May, J. 1996, A&AS, 115, 81 
Cabrit, S., Codella, C, Gueth, F, & Gusdorf, A. 2012, A&A, 548, L2 
Cabrit, S., Codella, C, Gueth, F, et al. 2007, A&A, 468, L29 
Caratti o Garatti, A., Froebrich, D., Eislotfel, J., Giannini, T., & Nisini, B. 2008, 

A&A, 485, 137 
Caselli, R, Hartquist, T. W., & Havnes, O. 1997, A&A, 322, 296 
Cesaroni, R., Felli, M., Jenness, T., et al. 1999, A&A, 345, 949 
Cesaroni, R., Neri, R., Olmi, L., et al. 2005, A&A, 434, 1039 
Codella, C, Bachiller, R., & Reipurth, B. 1999, A&A, 343, 585 



10 



S. Leurini et al.: Evidence of a SiO collimated outflow from a massive YSO in IRAS 17233-3606 

Codella, C, Beltran, M. T., Cesaroni, R., et al. 2013, A&A, 550, A81 

Codella, C, Cabrit, S., Gueth, R, et al. 2007, A&A, 462, L53 

Dayou, F. & Balanfa, C. 2006, A&A, 459, 297 

Flower, D. R., Le Bourlot, J., Pineau des Forets, G., & Cabrit, S. 2003, MNRAS, 

341,70 
Flower, D. R. & Pineau des Forets, G. 2003, MNRAS, 343, 390 
Guillet, v., Jones, A. R, & Pineau des Forets, G. 2009, A&A, 497, 145 
Guillet, v., Pineau des Forets, G., & Jones, A. R 2007, A&A, 476, 263 
Guillet, v., Pineau Des Forets, G., & Jones, A. R 2011, A&A, 527, A123 
Gusdorf, A., Cabrit, S., Flower, D. R., & Pineau des Forets, G. 2008a, A&A, 

482, 809 
Gusdorf, A., Giannini, T., Flower, D. R., et al. 201 1, A&A, 532, A53 
Gusdorf, A., Pineau des Forets, G., Cabrit, S., & Flower, D. R. 2008b, A&A, 

490, 695 
Hirano, N., Liu, S.-Y, Shang, H., et al. 2006, ApJ, 636, L141 
Hollenbach, D. 1997, in lAU Symposium, Vol. 182, Herbig-Haro Flows and the 

Birth of Stars, ed. B. Reipurth & C. Bertout, 181-198 
Hunter, T. R., Testi, L., Zhang, Q., & Sridharan, T. K. 1999, AJ, 118, 477 
Klein, B., Philipp, S. D., Kramer, I., et al. 2006, A&A, 454, L29 
Lada, C. J. & Fich, M. 1996, ApJ, 459, 638 
Le Bouriot, J., Pineau des Forets, G., Flower D. R., & Cabrit, S. 2002, MNRAS, 

332, 985 
Le Picard, S. D., Canosa, A., Pineau des Forets, G., Rebrion-Rowe, C, & Rowe, 

B.R. 2001, A&A, 372, 1064 
Lee, C.-R, Ho, P T. R, Palau, A., et al. 2007, ApJ, 670, 1 188 
Lefloch, B., Castets, A., Cemicharo, J., & Loinard, L. 1998, ApJ, 504, L109 
Leurini, S., Codella, C, Zapata, L., et al. 2011a, A&A, 530, A12 
Leurini, S., Codella, C, Zapata, L. A., et al. 2009, A&A, 507, 1443 
Leurini, S., Hieret, C, Thorwirth, S., et al. 2008, A&A, 485, 167 
Leurini, S., Pillai, T., Stanke, T., et al. 2011b, A&A, 533, A85 
Lopez-Sepulcre, A., Codella, C, Cesaroni, R., Marcelino, N., & Walmsley, 

CM. 2009, A&A, 499, 811 
May, P W., Pineau des Forets, G., Flower, D. R., et al. 2000, MNRAS, 318, 809 
McKee, C. F & Tan, J. C. 2002, Nature, 416, 59 
McKee, C. F & Tan, J. C. 2003, ApJ, 585, 850 
McMullin, J. P, Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, 

in Astronomical Society of the Pacific Conference Series, Vol. 376, 

Astronomical Data Analysis Software and Systems XVI, ed. R. A. Shaw, 

F Hill, & D. J. Bell, 127 
Nisini, B., Codella, C, Giannini, T., et al. 2007, A&A, 462, 163 
Panagia, N. 1973, AJ, 78, 929 

Qiu, K., Zhang, Q., Beuther, H., & Yang, J. 2007, ApJ, 654, 361 
Sault, R. J., Teuben, P J., & Wright, M. C. H. 1995, in ASP Conf Ser 77: 

Astronomical Data Analysis Software and Systems IV, 433 
Schilke, P, Walmsley, C. M., Pineau des Forets, G., & Flower, D. R. 1997, A&A, 

321, 293 
Schoier, F L., van der Tak, F. F S., van Dishoeck, E. F, & Black, J. H. 2005, 

A&A, 432, 369 
Su, Y.-N., Zhang, Q., & Lim, J. 2004, ApJ, 604, 258 
Tafalla, M., Santiago-Garcia, J., Hacar, A., & Bachiller, R. 2010, A&A, 522, 

A91 
Wang, K., Zhang, Q., Wu, Y, & Zhang, H. 201 1, ApJ, 735, 64 
Wrathmall, S. A., Gusdorf, A., & Flower, D. R. 2007, MNRAS, 382, 133 
Zapata, L. A., Ho, P T R, Schilke, P, et al. 2009, ApJ, 698, 1422 
Zapata, L. A., Leurini, S., Menten, K. M., et al. 2008, AJ, 136, 1455 
Zhang, Q., Sridharan, T. K., Hunter, T. R., et al. 2007, A&A, 470, 269 



11 



S. Leurini et al.: Evidence of a SiO collimated outflow from a massive YSO in IRAS 17233-3606, Online Material p 1 



O 
O 
O 
rsi 





1 1 1 1 1 1 1 1 1 [ [ [ 1 

H2 & Si0(5-4) 


-36°09'00" 


^^^ 




/iB^'V 


-36°09'20" 


KWLJ- 




Vail^ 


-36°09'40" 


- v^ — ^> 




■■< 




• J» 


-36°10'00" 





17''26""43^ 42^ 



41- 



a [2000] 



Fig.A.l. Integrated emission of the SiO(5^) line in the 
SMA+APEX combined data. The blue and red contours are the 
blue- and red-shifted emission (wbi = [-50, -20] km s"', from 
lOcr (1 Jy/beam) in steps of lOcr; v^i - [+10, +39] km s"', from 
lOcr (1 Jy/beam) in steps of lOcr). The grey image represents the 
H2 emission at 2.12/im. 



Appendix A: Joint deconvolution of the SIVIA and 
APEX data 

We performed a joint deconvolution of the SMA and APEX 
observations in the image plane using the task feather 
in the Commo n Astronomy Software Applications (CASA, 
iMcMuUui et a l.L 12007). A snapshot of the results of joint decon- 
volution of SMA and APEX data is shown in Fig. lA.ll While 
the large-scale details are filtered out by the interferometer, the 
combined image restores the large-scale information, i.e. (1) the 
combined image faithfully reproduces the large-scale (extended) 
structure, (2) the joint deconvolution image also conserved the 
total flux as observed with the single-dish. We tested this by 
smoothing the combined image to the resolution of the single 
dish beam. We also attempted to combine the data sets in the 
UV plane. However, the flux calibration factor (ratio of single- 
dish to interferometer flux in the overlapping UV space) was 
unexpectedly high. Although we could not ascertain the reason 
for the large scaling factor, the smaller area mapped with the 
single dish relative to SMA area might contribute to the failed 
combination. 

In Sect. 13.31 we demostrated that the SMA recovers the 
whole flux observed with APEX in the velocity ranges used in 
the analysis in Sect. |4] and therefore the combination of the two 
datasets was not needed for our purposes. 



